Bose-Einstein Correlations in e+e Annihilation 

and e+e- ^ W+W" 

W. J. Metzger 

■ Radboud University, Nijmegen, Netherlands 

o ; 

. Abstract. Results on Bose-Einstein correlations in e+e^ hadrons are reviewed. 
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O INTRODUCTION 

^ : 

I Bose-Einstein correlations (BEC) are those correlations which arise as a consequence of 

^ ■ Bose symmetry, which leads to an enhancement of the production of identical particles 

^ . close together in phase space. In this talk I will review some of the results on BEC from 

2 ! e+e^ interactions. 

Q>^ ■ To study correlations among q particles, we begin with the inclusive ^-particle density, 
O '• 

O' Pq{Pu-,Pq) = - J (1) 

X ; 

^ ■ where Oq is the inclusive ^-particle cross section. These densities are normalized such 

Dh! that 

^' ^Pl{p)dp = {n) (2) 

(3) 



J p2{PhP2)dpidp2 = {n{n-\)) 

J p2{puP2:P3)dpidp2dp3 = (n(n-l)(n-2)) , (4) 

etc., and are related to the factorial cumulants, C, by 

pi{pi) = Ci{pi) (5) 

P2{PI:P2) = Ci{pi)Ci{p2)+C2ipi,P2) (6) 
P3{Pl,P2,P3)) = Ci{pi)Ci{p2)Ci{p3) 

+ L Cl{pi)C2{p2:P3) 

3 perms 

+ C3{puP2,P3) ■ (7) 



Talk given at Fourth Workshop on Particle Correlations and Femtoscopy, Krommz, Czech Republic, 
August 15-17, 2005. A shortened version of this report will appear in the proceedings. 
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The ^-particle correlations are then measured by pq — nf=i Ci For q > 2, the 

correlations are a sum of "trivial", i.e., arising trivially from correlations of smaller q, 
and "genuine" correlations. For example, for ^ = 3 the trivial correlations are given by 
^3 perms^i(^i)^2(P2,/'3) and the genuine correlations by C3. 
It is convenient to normalize p and C: 

nLiPi(P') 

In order to study BEC, and not other correlations, one usually studies the ratio of the 
above-defined Rq to the Rq that one expects in the absence of BEC. This is equivalent to 
replacing the product of single-particle densities in dH) by the ^-particle density expected 
when BEC is absent, pog, giving: 

R, = -^. (10) 

POq 

This ratio is usually regarded as a function of Q, where = — (qm)^ with Mq the 
mass of the q particles and m the mass of each particle. If the particles have identical 
4-momenta. Q = 0. For 2 particles, is simply the 4-momentum difference. Thus, e.g., 
2-particle BEC are studied using 

^2(2) = -^^, (11) 

where the subscript q has been suppressed. 

It can be shown in a variety of ways that Rq is related to the spatial distribution of the 
particle production [1, 2]. For example, assuming incoherent particle production and a 
spatial source density of pion emitters, S{x), leads to 

R2(Q) = i + \G{q)\^ . (12) 

where G{Q) = J <ix:e'^^S{x) is the Fourier transform of S{x) . Assuming S{x) is a Gaussian 
with radius r results in 

i?2(2) = H-,-eV (13) 

Customarily, an additional parameter. A, is introduced in dTSt : 

R2iQ) = I + Xe-Q'''. (14) 

This parameter is meant to account for several effects: 

• partial coherence. Completely coherent particle production would imply A = 0. 

• multiple sources. 

• particle purity, e.g., experimental difficulty in distinguishing pions from kaons. 
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The assumption of a spherical (radius r) Gaussian distribution of particle emitters 
seems unlikely in e+e^ annihilation, where there is a definite jet structure. However, 
we must keep in mind that the BEC only occur among particles produced close to each 
other in phase space. In a two-jet event, a particle produced in the core of one jet would 
not be "close" to a particle produced in the core of the other jet. The volume in which 
hadronization occurs may thus be larger than a sphere of radius r and is not necessarily 
spherical. 

Furthermore, no time dependence has been considered, i.e., the source has been 
assumed to be static, which is certainly wrong. 

A number of other parametrizations have been considered in the literature. Never- 
theless, in spite of the above-noted limitations, this Gaussian parametrization (fT4ll is 
the one most frequently used by experimentalists. When this Gaussian parametrization 
does not fit well, an expansion about the Gaussian (Edgeworth expansion [3]) can be 
used instead. Keeping only the lowest-order non-Gaussian term, exp(— Q^r^) becomes 
exp(— 2^r^) ■ [l -I- jfH2,{Qr)~\ , where is the third-order Hermite polynomial. 

In the interest of comparison of as much data as possible, I shall only consider results 
using the Gaussian or Edgeworth parametrizations here. 



EXPERIMENTAL DIFFICULTIES 

There are a number of experimental problems which affect the results on BEC and their 
interpretation. 

Particle purity influences the value of A . In studying BEC of pions, it is often assumed 
that all particles are pions, which is not true. For example, in Z-boson decays approx- 
imately 15% of the particles are not pions. This lowers the observed value of A. The 
value of r for BEC involving a particle from a long-lived resonance and one produced 
directly, or involving particles from different long-lived resonances, will be larger than 
for two particles produced directly. Consequently, the resulting enhancement at small Q 
in RiiQ) will be narrower, possibly narrower than the experimental resolution, which 
in turn reduces the observed value of A. On the other hand, the effect of short-lived 
resonances is to increase the observed value of r, since it takes into account, in some 
average way, the distance the resonance travels before decaying. In Z-boson decays only 
about 16% of charged pions are produced directly, while 62% come from short-lived 
(r > 6.7 MeV) and 22% from long-lived resonances. Particles from weak decays are 
produced far away from the others, so that there are no BEC between pions from the 
weak decay and the rest. This results in a smaller value of A . About 20% of Z-bosons 
decay to bb. Not only the values of A and r may be influenced by resonances, but also 
the shape of R2. One can attempt to correct for all of these effects, typically by using a 
Monte Carlo model. However, the correctness of the model is open to question. 

A second problem is the choice of the so-called reference sample, i.e., the sample for 
which po is the density. Since it is impossible to turn off Bose statistics, this sample does 
not exist. It must be created artificially. Common choices are unlike-sign pairs, Monte 
Carlo models, and mixed events. When studying BEC in like-sign pion pairs one can use 
unlike-sign pion pairs to form the reference sample. The main problem is that the unlike- 
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FIGURE 1. (a) R2 (here called Q for Z-decay data and Monte Carlo; (b) the ratio of R2 of the data to 
that of the Monte Carlo. 

sign pairs have resonances which the like-sign pairs do not. This is illustrated in Fig.[T^, 
where R2 calculated using an unlike-sign pair reference sample is shown for Z-decay 
data and for Monte Carlo. Structure in R2 due to resonances in the reference sample is 
clearly seen. One might think of correcting for this using Monte Carlo, but this will only 
work if the resonances are well-described in the Monte Carlo. This is clearly not the case 
here, as is seen in Fig.[lj5. The solution adopted is to exclude the affected regions from 
the fit of (O to the data. These regions, indicated in Fig.[Tl cover large ranges in Q, and 
one must ask whether larger, or additional, excluded regions are needed. 

Using Monte Carlo as the reference sample avoids the different resonances of like- 
and unlike-sign pairs, but is still not completely free of resonance problems. Reflections 
of unlike-sign resonances and multi-body resonance decays, e.g., 77'(958) tz^k^t] — > 
TZ^ TZ^ K*^ , Still affect po- But these effects are small compared to the use of unlike- 
sign pairs as the reference sample. Another problem with using a Monte Carlo reference 
sample is the description of fragmentation. In a Monte Carlo model, fragmentation is 
controled by some parameters whose values are chosen to give a good description of the 
data. However, the data contain BEC. If the Monte Carlo does not include a description 
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TABLE 1. Results of fits H with and 
without FSI (Coulomb and S = nn 
phase shifts) to OPAL Z-decay data. 



FSI 


A 


r (fm) 


with 


1.04 ±0.03 


1.09 ±0.04 


without 


0.71 ±0.04 


1.34 ±0.04 



of BEC, the fragmentation parameters will be chosen so as to describe (in so far as is 
possible) BEC as well as the fragmentation itself. The reference sample will then have 
some enhancement in p at small Q, and its use as reference sample will result in a smaller 
value of A . The value of r may also be affected. On the other hand, a BEC model can 
be included in the Monte Carlo and the fragmentation parameters determined. Then the 
BEC model can be turned off to construct the reference sample. How good this reference 
sample is depends on how well the BEC were modeled. To model them well would 
require not only knowing how to correctly model BEC but also knowing the results of 
the BEC analysis before doing it. 

Event mixing involves taking pairs of particles from different events. For each real 
event, a similar event is constructed by replacing each particle by a particle from 
a different event. This procedure, of course, removes all correlations, not just BEC. 
Therefore, care should be taken in choosing the particles in order to have a reference 
sample similar to the real sample in characteristics such as jet structure. In this way some 
kinematic correlations can be preserved. Other correlations can be reinstated using the 
results of applying the same mixing procedure to Monte Carlo, i.e., by using the double 
ratio Ri/I^'^ instead of itself. Thus this method is also not entirely free of Monte 
Carlo uncertainties. For 2-jet events one can mix hemispheres instead of events. In this 
approach the reference sample is constructed by pairing a particle from one hemisphere 
with a particle from the other hemisphere which has been reflected through the origin. 

A third problem is final-state interactions, both Coulomb and strong. Identical charged 
particles will be repulsed by the Coulomb interaction, decreasing the observed value of A 
and increasing the observed value of r. The most common way to take this into account 
is to multiply Rq by a correction factor before fitting ([T4b to the data. The correction 
factor most commonly used is the so-called Gamow factor [4], which is the modulus 
square of the non-relativistic Coulomb wave function at the origin. For 2-particle BEC, 
given the resolution of the lep experiments, the correction is only a few per cent in the 
lowest Q bin (twice that if an unlike-sign reference sample is used), and is therefore 
frequently ignored. However for 3-particle BEC it is of the order of 10%. 

A different approach is to incorporate both the Coulomb and the 5 = TZK phase 
shifts into the wave function and derive a corrected formula for R2. An example of this 
approach is shown in Fig.|2l|@]. Significant differences are found in the values of X and r, 
which are listed in Table [l] However, such sophisticated approaches have not been used 
by any of the experimental groups. 

Finally, there is the effect of long-range correlations not adequately taken into account 
by the reference sample. As observed, e.g., in Fig. |2l Ri is not constant at large Q. 
To account for this in fitting the data, the right hand side of ([T4b is multiplied by an 
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FIGURE 2. Fits with and without FSI (Coulomb and 5 = KK phase shifts) to OPAL Z-decay data 01. 

appropriate factor, usually a linear dependence on Q, 

R2{Q) = r{i+^e-Q"'-")[\ + 5Q) , (15) 

although a quadratic dependence may also be used if necessary, and the normalization, 
7, is usually left as a free parameter. 

EXPERIMENTAL RESULTS 

Now let us turn to a comparison of results from various e+e^ experiments. 

2-particle BEC 

Dependence on the reference sample. The values of X and r found using identical 
charged-pion pairs from hadronic Z decays in the lep experiments, ALEPH |0, 0], 
DELPHI [8], l3 [9] and OPAL J 1, 12] are displayed in Fig.|3l 

Although these points are all supposed to be measurements of the same quantitities 
and are all determined by fitting the Gaussian parametrization to R2, there is little 
agreement among the points. Clearly there must be large systematic uncertainties not 
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FIGURE 3. A and r at = A^z found in the LEP experiments ||S H S 113 El El ■ Open arrow- 
heads at the end of error bars indicate that the uncertainties are statistical only, filled arrow-heads that the 
uncertainties are combined statistical and systematic. 



accounted for. Indeed, most of the points have only statistical error bars. However, some 
have error bars including also systematic uncertainties. But these are clearly insufficient. 

Solid points are corrected for pion purity, while open points are not. It is apparent that 
this correction increases the value of X but has little effect on the value of r. 

All of the results with r > 0. 1 fm were obtained using an unlike-sign reference sample, 
while all of those with smaller r were obtained with a mixed reference sample. It is clear 
that the choice of reference sample has a large effect on the values of the parameters. 
In comparing the results of different experiments, we must therefore be sure that the 
reference samples used are comparable. 

Dependence on the center-of-mass energy. The values of r found using identical 
charged-pion pairs in e+e^ annihilation is shown versus the center of mass energy, 
in Fig.m Keeping in mind that we should compare only results using the same reference 
sample, we conclude that there is no evidence for a ^/s dependence. 

Dependence on the particle mass. It has been suggested, on several grounds id, 
that r should depend on the mass of the particle as r oc 1/ ^Jm. Results from LEP 
experiments on r from 2-particle BEC for charged pions and for kaons, as well as the 
corresponding Fermi-Dirac correlations for protons and lambdas, are shown in Fig. |5l 
Again, restricting our comparison to results which have used the same type of reference 
sample (in this case mixed), we see no evidence for a 1/ ^/m dependence. Rather, the 
data suggest one value of r for mesons and a smaller value for baryons. The value for 
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baryons, about 0.1 fm, seems very small; if true it is telling us something unexpected 
about the mechanism of baryon production. 

Dependence on the particle multiplicity. The values of A and r from charged-pion 
2-particle BEC depend on the charged particle multiplicity, Hch' of the events. As seen 
in Fig.|6l A decreases with n^h while r increases. However, a similar dependence is also 
seen on the number of jets. When 2-jet events are selected, little, if any, dependence on 
Hch remains. For 3-jet events, r seems independent of n^h, although A does decrease with 
Rch- Thus, the dependence of A and r on n^h seems to be largely due to a dependence on 
the number of jets. 

Elongation of the source. The Gaussian parametrization (fT?b assumes a spherical 
source. Given the jet structure of e"'"e~ events, one might expect a more ellipsoidal 
shape. To investigate this, the parametrization is generalized to allow different radii 
for the direction along and perpendicular to the jet axis. The analysis is perfomed in 
the longitudinal center of mass system (LCMS), which is defined as follows: The pion 
pair is boosted along the jet-axis (taken, e.g., as the thrust axis), to a frame where the 
sum of the longitudinal momenta of the two pions is zero. The transverse axes, called 
"out" and "side" are defined such that the out direction is along the vector sum of the 
two momenta, pi+P2, and the side direction completes the Cartesian coordinate frame. 
This is illustrated in Fig. 13 
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FIGURE 5. Dependence of r on the mass of the particle as determined at -Js = Mz from 2-particle BEC 
for charged pions SHU 0113 [Hip, chai-ged kaonsjlljll and neutral kaons i3 El El and from 
Fermi-Dirac correlations for protons f2(f\ and lambdas Il22ll23ll . Open arrow-heads at the end of error bars 
indicate that the uncertainties are statistical only, filled arrow-heads that the uncertainties are combined 
statistical and systematic. The curves illustrate a 1 / y/m dependence. 
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FIGURE 6. Dependence of A and r from charged-pion 2-particle BEC on the charged particle multi- 
plicity, n^h, and on the number of jets 1.1 lil . 
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FIGURE 7. The definition of the longitudinal center of mass system (LCMS). 



In the LCMS, 




eL+eide+eout-(A£)2 

eL + esMe + eL(l-/3'), where /3 



Poutl + Pout 2 

£1+^2 



(16) 



The advantage of the LCMS is immediately apparent: The energy difference, and there- 
fore the difference in emission time of the pions, couples only to the out component, 
QoMt- Thus 2l and Qside reflect only spatial dimensions of the source, while gout re- 
flects a mixture of spatial and temporal dimensions. The Gaussian parametrization (fTSb 
becomes 



Assuming an azimuthally symmetric Gaussian shape, there is only one non-zero off- 
diagonal term, and G is given by 



out consistent with zero, and is then fixed to zero in the fits. Two-dimensional LCMS 
analyses have been performed by aleph |7] and DELPHI ll25ll . in which the out and side 
components are replaced by a transverse one, 2^ = -f Qgi^g. However, the inter- 
pretation of the corresponding parameter, rt, as a transverse radius is not unambiguous, 
since it includes the effect of the difference in time of emission. Both l3 and aleph fit 
not only (flTt . but a similar expression where G is replaced by a lowest-order Edgeworth 
expansion, which gives a better fit {e.g., in l3, a confidence level of 30% compared with 
3% for the Gaussian fit). The ratio of transverse to longitudinal radii found in the four 
experiments are shown in Table|2l The longitudal radius is clearly about 20% larger than 
the transverse radius. 

In Fig. [HI we see that the amount of this elongation increases when narrower 2-jet 
events are selected. 



i?2(eL,eside,eout) =7(1+ ^G)(1 + 5<2l + eeout + reside) • 



(17) 
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TABLE 2. Ratio of transverse or side radius to longitudinal radius from charged-pion 
2-particle BEC analyses of retail organisation 



Experiment 


Data 


Reference 


Gauss / 


2-D 


3-D 






sample 


Edgeworth 






DELPHI 


2-jet 


mixed 


Gauss 


0.62±0.02±0.05 




ALEPH 


2-jet 


mixed 


Gauss 


0.61±0.01±0.?? 






2-jet 


+ - 


Gauss 


0.91±0.02±0.?? 






2-jet 


mixed 


Edgeworth 


0.68±0.01±0.?? 






2-jet 


+ - 


Edgeworth 


0.84±0.02±0.?? 




OPAL 


2-jet 


+ - 


Gauss 




0.82±0.02±[J;0i 


l3 


all 


mixed 


Gauss 




0.80±0.02±[]?| 




all 


mixed 


Edgeworth 




0.81±0.02±[j:?^ 
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ycut 



FIGURE 8. The ratio of to rside for 2-jet events as function of the jet resolution parameter, jci 



It should also be mentioned that ZEUS ll26ll performed a similar 2-dimensional anal- 
ysis in deep inelastic ep interactions. The ratio rt/rL found, similar to that found by 
DELPHI and ALEPH, is independent of the virtuality of the exchanged photon. 

^r'^^r'^. In hadronization models with local charge conservation, e.g., string models, 
neutral pions can be produced closer together than identical charged pions. One could 
expect this to be reflected in a smaller BEC radius for n'^n'^ than for n^n^. Only two 
e+e^ experiments have attempted a BEC analysis of n^'^n^, l3 and OPAL f2§\. The 
experimental selections used in the two analyses are quite different, dictated as they are 
by the characteristics of the different detectors. While l3 requires the energy of the pions 
to be less than 6 GeV, opal demands the momenta to be greater than 1 GeV. Further, 
OPAL uses only 2-jet events, defined as having thrust larger than 0.9. 

The R2 distributions and fits are shown in Figs. |9l and [lOl respectively. In order to 
make a comparison with charged pions, l3 also analyzes n^n^ with a selection similar 
to its 7t^ selection. The resulting R2 is also shown in Fig. |2| The results of fits to these 
distributions, as well as two other n^n^ results are listed in Table|3l 

Comparison of the l3 values of r for n^n^ and n^n^ with the same selection indicates 
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FIGURE 10. R2 distribution for 7:°K° 



that the radius is smaller in the n n case, but the significance is only about 1 .5 standard 
deviations. A is also smaller in the n^n'^ case, with similar significance. 

For OPAL, the comparison is more difficult to make, since opal's charged-pi results 
use a different reference sample and different selection. Other experiments have found 
that the ratio of r using a mixed reference sample to that using unlike-sign is about 
0.68 (ALEPH l?]) or 0.56 (Delphi Applying such a factor would lower the opal 
value to about 0.62, which agrees well with the OPAL result for 71*^71*^. However, it is not 
clear what the effect of the 2-jet and ;r*^-momentum cuts is. In the l3 case, requiring the 
pions to have E < 6 GeV and using a Monte Carlo reference sample rather than a mixed 
one decreased r by about 30%. It is therefore conceivable that the OPAL requirement of 
p> 1 GeV would increase r, in which case r for Tt^Tt^ would be smaller than for k^k"^. 
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TABLE 3. Results of fits to R2 for n^n'^ and tt^tt* jHllIll. The indicated 

uncertainties combine statistcal and systematic uncertainties. 



Experiment 


Selection 


Reference 
sample 


r(fm) 


A 


n^n^ l3 

OPAL 


<6GeV 
p„ > 1 GeV, 2-jet 


MC 
mix 


0.31±0.10 
0.59 ±0.09 


0.16±0.09 
0.55±0.14 


n'^n'^ l3 
l3 

OPAL 


Ej,<6GeV 


MC 
mix 
+- 


0.46 ±0.01 

0.65 ±0.04 

1 00+^-^^ 
^■^^-0.10 


0.29 ±0.03 
0.45 ±0.07 
0.57 ±0.05 



3-particle BEC 

There have also been analyses of BEC among three charged pions. As mentioned 
in the Introduction, with three (or more) particles the correlations may be classified as 
trivial, simply a consequence of lower-order correlations, and genuine. We study these 
correlations using RsiQs) = piQs) /poiQs)- Note that Qj = + Q23 + Qli- The same 
assumptions that lead to (O for R2, lead to 1E9II30I1 

RliQij) = l + X\Giqij)\ (19) 
^3(212,223,231) = l+A(|G(2l2)|'+|G(223)P+|G(23l)P) 

. ' 

from 2-particle BEC 

+ 2Ai-59^{G(2i2)G(223)G(23i)} , (20) 

^ V ' 

from genuine 3-particle BEC 

where G{Q) = J dxe'^^S{x) = \G\ e"^ is the Fourier transform of the source density, and 
9? denotes the real part. Note that R3, unlike R2, depends on the phase of G{Q). We 
define R^^^ as the R3 that would occur if there were no 2-particle BEC: 

i?r(2i2, 223,231) = l+2Ai-59t{G(2i2)G(223)G(23i)} ■ (21) 

Defining 

CO = cos(<|)i2 + + (1>13) , (22) 

the above equations yield 

^r(23) - 1 

CO = — , , (23) 

2v/ (^2(212) - l)(i?2(223) - l)(i?2(2l3) - 1) 

which simplifies in the case of a Gaussian to 

^r(23)-i 

CO = ^, ^ ' . (24) 

2v/^2(23)-l 

If the particle production is completely incoherent, the phase (^tj is expected to be zero, 
and consequently we should find (O = 1. As we have seen, incoherence also implies 
A = 1, but A is affected by many other factors, which should not affect ft). 
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TABLE 4. Results of fits to R2 and Rf^" using Gaussian 
and Edgeworth parametrizations |9]. 



fit to 




Gaussian 


Edgeworth 


R2 
Rf 


X 


0.45 ±0.06 ±0.03 
0.47 ±0.07 ±0.03 


0.72 ±0.08 ±0.03 
0.75±0.10±0.03 


Ri 
Rf 


r 

(fm) 


0.65 ±0.03 ±0.03 
0.65 ±0.06 ±0.03 


0.74 ±0.06 ±0.02 
0.72 ±0.08 ±0.03 
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FIGURE 11. RiiQ) measured by l3 j^l- The full circles represent the data, while the open circles 
correspond to results from a Monte Carlo model without BEC. The dashed and full lines show the fits of 
a Gaussian and an Edgeworth expansion, respectively. 



The l3 measurements [9] of R2 and i?!^" in hadronic Z decays are shown in Figs.fTTl 
and[T21 The distributions are fit with both a Gaussian and an Edgeworth parametrization. 
The Edgeworth parametrization provides a better fit of the data. It is reassuring to note 
that the values of A and r obtained from the fits to R2 and R^^'^, which are listed in 
Table m agree perfectly. 

The dashed lines in Fig. [121 are the predictions of using (1211) with the results of 
the fits to R2 assuming (0 = 1. These dashed lines are quite close to the solid lines, which 
represent fits to Rf^^iQs), confirming that (O = 1 is a reasonable hypothesis. 

Another way to examine this hypothesis is to compute ft) using (l23t for each bin in 
Qj, using the measured ^^^^"(Qs) and R2 from its fits. The results are shown in Fig.fTBl 
We conclude that O) is perfectly consistent with unity, and consequently that particle 
production is completely incoherent. 

BEC among three pions have previously been observed at lower energies iflsL [l4ll. 
At LEP, DELPHI iHl] and OPAL Jsl] have also observed genuine 3-pion BEC. With the 
exception of DELPHI, which did not do a comparable 2-pi analysis, they all report that 
the values of A and r obtained for 2- and 3-pions are consistent. Unfortunately, none of 
these experiments performed an analysis using O). 
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FIGURE 12. i?3'^"(23) measured by l3 101 ■ The full circles represent the data, while the open circles 
correspond to results from a Monte Carlo model without BEC. The full line in (a) is a fit of a Gaussian 
parametrization, while in (b) it is a fit of an Edgeworth parametrization. The dashed lines show the 
expectation from i2l\ with the results of the fits to R2 assuming (O — I. 

BEC in hadronic W decays 

W qq. Having found no evidence for a center-of-mass energy dependence of 
BEC, and noting that in any case the masses of the Z and W are not much different, 
there is only one reason to expect BEC to be different in W decays than in Z decays, 
namely the different flavor composition. Whereas about 20% of Z decays are to bb, 
almost no W decays involve a b-quark. The long lifetime of the b results in diminished 
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FIGURE 13. The ratio co i23\ as function of computed using the measured /^^"(Sb) and ^2 from 
its fits. 



BEC. Thus, we should expect BEC in W-decay to be like that in the decay of the Z to 
light (udsc) quarks. This is indeed found to be the case ll33ll34ll . as is shown in Fig.fT?! 
where RiiQ) for hadronic decays of the W-boson produced in e+e~— >^ W+W^ — > qq£v 
is compared to that of Z-boson decays to all flavors and to udsc flavors only. 

Inter-string BEC. A more interesting case is e+e^^ W+W^ qqqq. The W- 
bosons are produced not far above threshold and consequently travel only about 0.7 fm 
before decaying. This is smaller than the distance over which hadronization occurs. 
Therefore one expects a significant degree of overlap of the two hadronizing systems, 
resulting in BEC not only between particles produced by the same W, but also between 
particles produced by different W-bosons. On the other hand, in the string picture no 
BEC is expected between particles arising from different strings. 

Whether or not this so-called inter-string BEC exists is thus a fundamental question 
for the string picture. It is also important for the measurement of properties of the W- 



boson, in particular its mass 11351] . Improper simulation of inter-string BEC in IVlonte 
Carlo programs would lead to a bias in the mass measurement. 

All four LEP experiments have studied this question [33, 36, 37, 3^ 34, 3^. The basic 



method [40] is to test the expectation of no inter-string BEC. If the two W-bosons decay 
completely independently, the two-particle density in the 4-quark channel is given by 

pAq{p\,P2)= p+iPhPi) 1,2 from W+ 

+ p^{pup2) 1,2 from W" 

+ p+{pi)p-{p2) 1 fromW+, 2fromW^ 

+ p+{p2)p^{pi) 1 fromW", 2fromW+. 
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FIGURE 14. R2iQ) measured by l3 fsT] for hadronic W-boson decays from e+e 
and hadronic Z-boson decays to all flavors and to udsc flavors only. 



qqiv 



Assuming p+ = p = p2q, which would only be strictly true if there is complete overlap, 

p4q{Pl,P2) = 2p2q{pi,P2) + 2p2q 1 ) p2q (/'2 ) • (25) 



The density p4q(pi,P2) is measured in e+e^^ W+W^ qqqq and p2q(pi,P2) 
in e+e^^ W+W^ qqiv. The remaining term, p2q{pi)p2q{P2), is estimated by 
Prmx{p\TP2) obtained by mixing £+vqq and qq£^v events after removal of the 
and £~ . Thus, (l25t becomes 



p4q(e) = 2p2q(e)+2pmix(e) ■ 

Various quantities are defined to test the validity of ( E^ : 

Ap(2) = p4q(e)- [2p2q(;?l,P2)+2pmix(;?l,P2)] 

m) = 



p4q{Q) 



2p2q(e)+2pmix(e) 
Ap(e) 



(26) 

(27) 
(28) 

(29) 



2pmix(e) ■ 

The quantity 5i{Q) is actually the correlation function of genuine inter-W BEC ll4lll and 
is thus not only a test of no inter-W BEC, but the quantity of interest if inter-W BEC 
indeed exists. 
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FIGURE 15. di{Q) measured by Delphi js^ for (a) like- and (b) unlike-sign pion pairs. Also shown 
are the expectations of Monte Carlo models including BEC among all pions or including BEC only among 
pions from the same W. 



All four LEP experiments have now reported final results 3^ 34, [39']. Fig. [T51 
shows results of Delphi 113 8h . the only experiment which claims to have seen signif- 
icant inter-W BEC. The 5i distribution of like-sign pions (Fig. [T3b ) clearly shows an 
enhancement at small Q, while no enhancement is seen in the Monte Carlo distribution 
where BEC, as modeled by the BE32 model f4^, is included only between pions from 
the same W. When pions from different W-bosons are modeled with the same BEC as 
pions from the same W, a larger enhancement is seen than that of the data. 

However, the interpretation of the enhancement in the data as inter-W BEC is some- 
what clouded by the observance of an enhancement, albeit smaller, in the 5i distribution 
of unlike-sign pions (Fig.[T3b). In this distribution the enhancement in the data is larger 
than that for Monte Carlo with inter-W BEC. 

The final results for the four LEP experiments are compared in Fig. [T^ Here the 
results are expressed as the ratio of the effect seen to that expected in the BE32 model 
with the same inter-W BEC as intra-W BEC. The measurements indicated by arrows 
are combined to give the preliminary LEP result 11431] of 0.17 ±0.13. The data are thus 
compatible with no inter-W BEC. 



BEC in the Lund string model 

In the Lund string model, the longitudinal break-up of the color string is governed 
by the area law. The matrix element to get a final state depends on the area A: M = 
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FIGURE 16. The amount of inter-W BEC observed expressed as a fraction of the prediction of the BE32 
modeL The measurements indicated by arrows are combined to give the preliminary LEP result i43ll . 



exp[{iK — b/2)A], where K is the string tension and b is a decay constant with values 
K" ^ 1 GeV/fm and ~ 0.3 GeV/fm. Consider the string break-up illustrated by the solid 
line in Fig.[T7J which spans the area A. Suppose that the particles 1 and 2 are identical. 
The final state produced by interchanging them would be produced by the string break- 
up with area A + AA. Transverse momentum arises via a tunneling mechanism, which is 
also related to Z?. To incorporate BEC in the string model, the probability of a final state 
should be taken as the square of the sum of the matrix elements corresponding to the 
areas of all the permutations of identical bosons |44, 45, 46]. This model results in BEC, 
including genuine 3-particle BEC. It predicts that the longitudinal radius is greater than 
the transverse radius and that the radius for neutral pions is smaller than that for charged 
pions. In such a model there is no mechanism for creating correlations between particles 
from different strings. 

The model has been incorporated in Monte Carlo for a qq string. While the formalism 
to do so for the more realistic case of a string with multiple gluons has been worked out 
i47.] . a successful Monte Carlo implementation has, unfortunately, so far proved elusive. 



CONCLUSIONS 



We have seen that the study of BEC in e+e presents a number of problems, both ex- 
perimental and theoretical. Consequently, values obtained for parameters vary consid- 
erably among experiments, even when the same parametrization is used. Nevertheless, 
certain features are clear: BEC, both 2-particle and genuine 3-particle, exist; they seem 
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FIGURE 17. The break-up of a Lund string into hadrons, showing two break-up patterns giving the 
same final state with two identical particles interchanged.. 



independent of center-of-mass energy; the source shape is somewhat elongated in the jet 
direction; and the (Fermi-Dirac) radius for baryons is smaller than the radius for mesons. 
Experimentally, it is not clear whether the radius for neutral pions is smaller than that 
for charged pions. BEC in W decay is the same as in light-quark Z decay. The data are 
compatible with no inter- W BEC. 

The implementation of BEC in the Lund string model appears consistent with these 
experimental findings. However, the experimental evidence that pion production is com- 
pletely incoherent seems at odds with the coherent addition of amplitudes in the model. 
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